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Abstract 

We show how the separability problem is dual to that of decom- 
^ , posing any given matrix into a conic combination of rank-one 

partial isometries, thus offering a duality approach different to 
the positive maps characterization problem. Several inmediate 

■ consequences are analyzed: (i) a sufficient criterion for separabil- 
, ity for bipartite quantum systems, (ii) a complete solution to the 

lO ■ separability problem for pure states also of bipartite systems indc- 

pendent of the classical Schmidt decomposition method and (iii) 
^ c| a natural generalization of these results to multipartite systems. 
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C7 1 1 Introduction 



Entanglement is at the heart of the most surprising quantum phenomena, es- 
pecially those sustaining the quantum processing of information (cf. e.g. P3G]). 
However, a complete characterization of entanglement of compound quantum 
systems is not nowadays at hand. There exists a great deal of either sufficient 
or necessary conditions for a system to show entanglement (cf. 0J [SJ E3 E| 
and multiple references therein), and those both necessary and sufficient criteria 
either apply to a limited number of systems [HIE] or lack of an operational (user- 
friendly) character |1U1 111! 1121 . Among all of them, the most celebrated is 
the well-known PPT or Peres-Horodecki criterion [OHEI; which discerns between 
entangled and separable states just by checking the positivity of the density ma- 
trix of the bipartite compound system after performing the partial transposition 
1 . The grandeur of the PPT criterion resides on the translation of the separa- 

1 Given a matrix Q £ JV[ m (g) M n , which can always be written as Q = 53ij=i ^ij ® Qij, 
where Eij denotes the Weyl matrices, the partial transposition with respect to the first factor is 
defined as Q Tl = J3y = i ^Tj^Qij an d with respect to the second factor as Q T2 = J3y = i ^ij® 
Qij, where A T , as usual, denotes the transpose of A. 



1 



bility question to a dual problem, namely, the characterization of positive linear 
maps between matrix algebras jHj . Since this long-standing mathematical prob- 
lem has only been solved for maps between two two-dimensional |15j and a 
two-dimensional and a three-dimensional matrix algebra , this criterion only 
works for 2x2 and 2x3 systems. The ultimate reason for this limitation then 
turns out to be the lack of knowledge about the structure of the set of positive 
linear maps. In the multipartite case, the situation is even more dramatic, since 
no full characterization is known. 

Here we formulate an equivalent dual problem to that of separability of quan- 
tum states by resorting to two mathematical tools commonly used in Matrix 
Analysis (cf. e.g. ^Tj), namely the vec operator carrying matrices onto vectors 
and the tilde transform, realigning matrix elements into another matrix with the 
same total number of elements, although in general, different dimensions. As a 
direct consequence we obtain both a sufficient condition for separability and a 
characterization of product states for bipartite and multipartite systems of any 
dimension, thus offering an alternative to the classical Schmidt decomposition 
in the case of pure states. Very close objects have already been used in this con- 
text in the so-called realignment method |18j . although the use of a (Ky-Fan) 
norm in this approach drives us only to a necessary condition for separability 
(cf. ^3 |201 for an equivalent approach and [2] for the generalization of this 
method to multipartite systems.). 

We report our results as follows. In sectionElwe introduce the mathematical 
tools, namely the vec operator and the tilde transform; in section|3we establish 
the new duality problem and analyse some of its first consequences for bipartite 
systems of any dimensions; in section 0] we extend preceding results to the 
multipartirte case; we illustrate the results with some examples in section [5] we 
conclude in section [5] with some comments and remarks. 

2 The vec operator and the tilde transform 

Throughout the whole paper we will denote the set of n-dimensional matrices 
by M n and vectors of any dimension by the usual bra-ket notation. The vec 
operator is commonly defined as a linear isomorphism carrying matrices M of 
arbitrary dimensions raxnto mn— dimensional vectors obtained by stacking 
the matrix columns successively from left to right (cf. [22] for a mathematical 
review) : 

Definition 1. Let Q = [q± . . . q n ] £ Ai n , where qi demotes the ith column of Q. 
Then the vec operator is defined as 
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vec : M n 



\vec{Q)) 



( 9i 



\ 9n 



(2.1) 



This stacking is conventional and generalizations of this operator can be 
found in the literature |23| with different purposes and applications. The gen- 
eralization to rectangular matrices is elementary, although useless for our im- 
mediate purposes. It is straightforward to prove the following 

Proposition 1. The vec operator is a linear isomorphism. 

Related to it, one can also introduce a so-called tilde transform [53], which 
for the same reasons as before we will only define for square matrices (the 
generalization to rectangular matrices is also elementary): 



Definition 2 

tilde transform of Q is defined as 



. Let Q € M m (£> M n , which can always be written as Q = 
ij, where Eij are the m— dimensional square Weyl matrices. The 



Q = t(Q) = 



\vecQ l] )(vecE\ j \ 



(2.2) 



Notice that Q is an n 2 x m? rectangular matrix. This definition drives us to 
a realigned matrix different from that used in |18j to establish a new (necessary) 
separability criterion, namely, the realignment criterion. It is straightforward to 
convince oneself that the realignment method can also be established with the 
above convention, since the Ky-Fan norm properties |17j of the tz/rfe-transformed 
matrices are equivalent. Let us provide the following illustrative example in 
Mi <g> Mi: 



Q 



( 9ii 912 

921 922 

931 932 

\941 942 



913 914^ 

923 924 

933 934 

943 944 / 



t(Q) 



(q\\ 931 913 933^ 

921 941 923 943 

912 932 914 934 

\922 942 924 944/ 



(2.3) 



The main property why this definition has been adopted is the following 
Proposition 2. Let Q\®Qi G M m ®M n - Then t{Qi®Q 2 ) = \vecQ 2 )(vecQ\\. 
Proof. The proof follows inmcdiatcly from the definition of the tilde transform 
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and the linear properties of the vec operator: 
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t{Qi®Q 2 ) 



(2.4a) 



ij=i 



(2.4b) 



| W ecQ 2 )(«ec(£ gg 5 *^)! 



(2.4c) 



ij=l 



\vecQ 2 ){vecQ\\ 



(2.4d) 



□ 



As before, it is elementary to prove the following 

Proposition 3. The tilde transform is a linear isomorphism. 

Notice that in Quantum Mechanics we are mostly involved with Hermitian 
matrices, thus for any product state p\ ® p 2 , the latter proposition reduces to 
t{p\ ® p 2 ) = \vecp 2 }(vecpi\. Had we used the matrix realignment criterion of 
|18| . we would have ended with the property t(A <g> E) = \vecA)(vecB*\, which 
introduces the nuisance of taking complex conjugates of all elements of B. In 
table ^ we include all other possible definitions of a tilde transform driving to a 
similar tensor product property. We have chosen definition 8 in order to have 
either vecEij in the bra part of its definition (so that the rcalingment does not 
involve also complex conjugation) and the adjoint of one of the factors in the 
tensor product property (since this operation will not affect Hermitian matri- 
ces and, in particular, density matrices). Any of these definitions allows us to 
reproduce the subsequent results in a similar fashion. 

The matrices of the form \u)(v\ belong to a special class of matrices called 
partial isometries, thus we will refer to them as rank-one partial isometries |17j . 
Also, we will give vec -1 , the inverse of the vec operator, carrying vectors to 
matrices, a special name: mat = vec -1 . Its action is clear from that of vec 
(read relation l|2.1J) from right to left). Note that applying the tilde transform 
twice to a matrix leaves it invariant, i.e. the tilde transform is its own inverse 
(check for instance equation (|2.3(l ). 

Finally, let us comment that under mathematical rigor we should have de- 
noted the vec operator as vec n and the mat operator as mat n , keeping clear 
upon which matrix algebra they are operating. This is also valid for the tilde 
transform, which should have been denoted as t TO ,„. However in order to ease 
the notation, we will assume that these matrix vector spaces are fixed from the 
beginning so that the corresponding dimensions are known. 
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Definitions of t 



2- t(Q) = 



\vecEij){vecQij 
vecEij){vecQ\- 
vecEji) (vecQij 



Em 
ij — 1 

Em 
ij — 1 

Em 
ij — 1 

Em 
ij=l 

E 

Em 
ii=l 



m 

y=l 



vecQ*j) (vecEij 
vecQij)(vecEij 

■Qij)( vecE ij 



vecC 



vecQij)(vecEj 



Tensor product property 



t(Qi®Q 2 
t(Qi <8> Q 2 
t(Qi ® 
t(Qi ® 

i(Ql ® <2: 

t(Qi ® Qa) 
t(Q 1 ®Q 2 ) 
t(Qi 



) = 
) = 



2 = 



) = 



vecQi)(vecQ2 
vecQ{){vecQ 2 
vecQi) (vecQl 



vecQ\)(vecQ 2 
\vecQ* 2 ){vecQ\ 
vecQ 2 )(vecQl 
vecQ\) (vecQi 
vecQ 2 )(vecQ\ 



Table 1 : Definitions of the tilde transform with the corresponding tensor product 
property. 



3 The dual formulation and some consequences 

We formulate our main result for bipartite systems, which establish the dual 
formulation of the separability problem: 

Theorem 1. Let p be a density matrix of an m x n bipartite system. Then p is 
separable if, and only if, p admits a conic decomposition into rank- one partial 
isometries, i.e. 

K 

p = y^ j Hi\u i ){v i \ p t >0 (3.5) 

i=l 

where mat\ui) > and mat\vi) > for all i = 1, . . . , K . 

Proof. If p is separable, then there exist positive matrices pi G M. m and c, G 
M. n (density matrices, indeed) and positive numbers pi such that p = Ei=i P^iPi® 
o~i. Applying the tilde transform upon p, we obtain 

K 

p = Pit(pi <E> crj) = (3.6a) 

= jUj |«ecg-j) («ec/?t | = (3.6b) 
i=i 

= Hi\vecx7i){vecpi\ = (3.6c) 
i=i 
if 

= (3.6d) 

?:=i 

where clearly mat\ui) — cr, > and mat\vi) = pi > 0. 
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Conversely, if such a decomposition exists, then 

P = KP) = P = ( 3 - 7a ) 

K 

= ^int{\ui)(vi\)= (3.7b) 

K 

= fii(mat\vi})^ (g) mat\ui) = (3.7c) 

i=l 
K 

= Himat\vi) ® mat\ui) (3.7d) 

i=l 

Thus, since mat\ui) and mat\vi) are (Hermitian) positive matrices, p is separa- 
ble. □ 

As a first consequence of this result, the singular value decomposition J7j 
provides a method to find a conic decomposition into partial isometries, thus 
furnishing a sufficient criterion of separability: 

Corollary 1. Let p be a density matrix of an m x n bipartite system. Let 
{\vi)}i=i,..., q and {\wi)}i = i.,,,^ q be left and right singular vectors, respectively, 
associated to non-null singular values of p. If mat\vi) > and mat\wi) > for 
all i = 1, . . . , q, then p is separable. 

Proof. Let us recall |17j that the singular value decomposition of a matrix Q = 
VY,W^ can be rewritten as Q = <J i\ v i){ w i\: where \vi) and \wi) are left 

and right singular vectors associated to non-null singular values o~i. The rest 
follows elementary from theorem ^ □ 

The necessity is unattainable in general, since singular vectors are always 
orthonormal, i.e. (vi\vj) = (wi\wj) = <Sy. This restriction is, however, absent 
from theorem ^ thus the singular value decomposition cannot provide the re- 
quired generality. Nevertheless, a notable advantage of this sufficient criterion 
is that, in the positive case, it yields by construction a convex combination of p 
in product states. 

Another inmediate consequence of our result is a necessary and sufficient 
criterion to know whether a given state is a product state or not: 

Corollary 2. Let p be a density matrix of a bipartite system. Then p ~ p\® pi 
if, and only if, rank(p) = 1 and the unique left and right singular vectors \u) and 
\v) associated to the non-null singular value of p satisfy mat\u) > 0, mat\v) > 0. 

Notice that this result offers a separability check for pure states alternative 
to the common Schmidt decomposition |25l I26j , which, in addition, in the pos- 
itive case also provides the corresponding factors. Starting from a pure state 
|^}, all we have to do is to apply corollary [3 to p^ = |^')( 1 I r |, together with the 
further conditions r(mat\u)) = r(mat\v)) = 1. Moreover, as shown below, this 
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method can also be exported to the multipartite case in a natural way. 

Another characterization of those relevant rank-one partial isometrics fulfill- 
ing the conditions of theorem ^ can be obtained by resorting to the structure 
of the cone of positive matrices Q > if, and only if, there exist positive 
numbers and orthogonal projectors P Ui = \ui)(ui\ such that Q = AjP Ui . 
Keeping in mind that vecP u = \u*) ® \u), then theorem ^ can be reformulated 
as 

Theorem 2. Let p be a density matrix of an m x n bipartite system. Then p is 
separable if, and only if, p admits a conic decomposition into rank- one partial 
isometries, i.e. 

K 

p = ^2ni\ui}(vi\ pt>0 (3.8) 
i=i 

where for each i = 1, . . . , K 

\m) = \x*) ® \xi) (3.9a) 

h) = \yt)®\yi) (3.9b) 

Proof. Starting from p = J^. \i\ui){vi\, where Xi > and mat\ui), mat\vi) > 0, 
insert \ Ui ) = Ek^l^*) ® 14^) an d |«i> = E^fl^*) ® 1^) to obtain 

p = J2^^ (l4 l) *)»l4 4) )) ((^1® (V?\) (3-10) 

Merge the enumerable indices into a single one to arrive at the desired result. 

□ 

Regretfully this result does not either provide us with an operational separa- 
bility check, since an algorithm to find such a decomposition into the preceding 
rank-one partial isometrics is in general unknown. However, a slight difference 
with respect to the positive maps characterization problem is detected. Let us 
make use of the language of the theory of cones [2U 2 ■ The positive maps char- 
acterization problem can be formulated in terms of the analysis of the structure 
of the cone IT m „ of positive maps between an m- and an n-dimensional ma- 
trix algebra mapping positive matrices of M. m into positive matrices of A4 n : 
the set of extremal elements of n TOi „, denoted as Extn mi „, is only known in 
the cases (m,n) = (2,2), (2,3) and (3,2). The reformulation proposed here 
reduces the separability problem to find a criterion to discern whether a given 
arbitrary n 2 x m 2 rectangular matrix belongs to the cone K n 2 m 2 1 whose set of 

2 A cone K is a closed set K of a vector space V such that \K C K for all A > 0; in the 
theory of cones referred to above a cone is always (i) convex (\K + (1 — X)K C K), (ii) pointed 
K fl (-K) = {0} and (iii) reproducing K — K = V . One of the most famous cones is that 
of positive semidefinite matrices (cf. I27I V The set of cone-preserving linear maps can be also 
proved to be a cone. Among the most relevant properties of a cone is that any of its element 
can be expressed as a linear combination of its extremal elements with positive coefficients 
(conic combination). 
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extremal elements is, by definition, Exti^„2 m 2 = : \u) = \x*} <£> \x) £ 

C n ® C", \v) = \y*) ®\y) £ C m <g> C m }. Now the set of extremal elements is 
known, thus everything is reduced to find an operational criterion to decide 
whether p belongs to K n i m i or not. 



4 The multipartite case 

To analyse the multipartite case, we need to generalize the tilde transform. Let 
us first recall that an arbitrary matrix Q £ M. ni SS • • • ® -M.n m can be written 
in any of these m forms: 

712 n m 

Q = E E Q£L.A mjm ® E hh®---®E imjm (4.11a) 

*2 ,5*2=1 im,im=l 

= E E E E nn®Qili 3j3 ...i mjm ®---®E imjm (A.llb) 

il,il = l i3,i3 = l Cjm^l 

rii n m — 1 

»l,Jl = l in TO -l,Jn TO -l=l 

(4.11c) 

where j denotes the Weyl matrices in the space M np - The matrices Q^ k > = 
Q^jj have dimension xji^. A systematic procedure to find these matrices 
Q( k > arises from the use of the vec-permutation matrix P(m, n) = Y^T=i £j=i E ij® 
Eji and the property B Cg> A = P T (m, n)(A ® B)P(m, n), where A £ A4 m and 

b 6 M n 113. 

With the same spirit as before we introduce the generalized tilde transforms: 

Definition 3. Let Q £ A4 ni 8> ••■ ® -Mn m - ?7ie /ci/i iiZcfe transform of Q is 
defined as 

ni nk—i n k+i n m 

w)=E- E E ■•• E i^..., m JH< J1 ®--'<_ 1Jt . 1 < +1]t+1 ®-< J j 

(4.12) 

This generalized tiWe transform also satisfies a similar property: 
Proposition 4. Let Q £ M ni ® ■ • • ® M nm . Then 

m 

*fc(0<9») = |«ecQ fc )(«ec(Qt®---®Q^_ 1 «>Q^ +1 <8>---ig)Q^| (4.13) 

2 = 1 
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Proof. It is an inmediate generalization of the proof of proposition [21 □ 

The analysis of the separability properties of a multipartite state depends 
on the particular partition chosen for the compound system. Here we will focus 
on full separability i.e. on the conditions upon which a state p of an m— partite 
system can be written as p = i Ai 1 ...i m pi 1 ® ■ ■ ■ ® p% m . 

Theorem 3. Let p be a state of an ri\ x • • • x n m m— partite system. Then p is 
fully separable if and only if tk(p) admits a conic decomposition into rank-one 
partial isometries 

qk 

t k (p) = J2^\u^)(v^\ (4.14) 

i=l 

where mat\u^) > and mat\v^) > for all i = 1, . . . ,qi~ and for all k = 
l,...,m. 

The proof is a tedious, although elementary, algebraic generalization of that 
of theorem ^ The preceding corolaries can also be adapted to the multipartite 
case, in particular, the use of the singular value decomposition also allows us to 
define the following sufficient criterion of full separability: 

Corollary 3. Let p be a state of an ni x • • • X n m m— partite system. Let 
,qk an d i\ w i ! )}i=i, ^ e l e ft an d right singular vectors, respec- 
tively, associated to non-null singular values of tk(p). // mat\v^) > and 
mat\w^) > for all i = 1, . . . , qk and for each k = 1, . . . , m, then p is fully 
separable. 

Once more, the orthonormality of singular vectors prevents us from attaining 
the necessity of this criterion. As above, we can also establish a necessary and 
sufficient test of separability for pure states, thus circunvemting the lack of 
generalization of the Schmidt decomposition to multipartite systems |28) . 

Corollary 4. Let |^} be a pure state of an n\ X • • • X n m m— partite system. 
Then \*S>) is a product state (fully separable) if and only if rank(tk(p\s/)) = 1 for 
all k = 1, ... ,m and the unique left and right singular vectors |t>W) and |u/ fc '} 
corresponding to the non-null singular value of tk(py) satisfy mat\v^ k ') > 0. 
matlw^) > and r(mat\v^)) = r(mat\w^ k '}) = 1 for each k = 1, . . . , m. 

The analysis of other partitions can be easily undertaken with the following 
recipe. For concreteness' sake, let us suppose we want to investigate whether 
a given 5— partite system pure state \$>) can be decomposed as |\&) = |V'}i2 ® 
|0}3 ® W)a5, i-e. whether it admits a (12)(3)(45) separable partition. Then we 
treat systems 1 and 2 as an n\ni single system, as well as systems 4 and 5 as 
an 7x4715 single system, and then apply the preceding result. The analysis of the 
separability of the factors \tfj) 12 and |"0}45 can also be undertaken with the same 
tools, thus providing a systematic method to detect multiseparability. These 
results complement those of [55] EDI ■ 



5 Examples 



As a first elementary and well-known example we will prove how the completely 
depolarized state in an m x n system is separable: 



1 1 m 1 

P = Iran = V E^> g <%I„ p = \vecl n ) (vecl m \ (5.15) 

mn mn ^— < y J mn 

*j 

Since mat\vec\k) > for all k = 1, 2, . . . , p is separable. 

A second example is provided with a Bell basis state, which serves to illus- 
trate that if a state is entangled, then it cannot fulfill the conditions of corollary 
^ We firstly need the following elementary property of mat: 

2 

Proposition 5. The mat operator is an isometry between C l and A4 n , where 
the inner products are the standard complex scalar product (u\v) = J2k u k v k an< ^ 
the trace scalar product (A\B) = tr(A^ B). 

Then, it is easy to convince oneself that the singular value of any Bell basis 
state is 1 (fourfold) . and since it is impossible to find four pairwisc orthonormal 
vectors vi € C 4 such that matvi > and V = [i>i V2 V3 V4] is unitary 3 , we can 
conclude that the conditions of corollary 2] are impossible to fulfill. 

Another example illustrating how the singular value decomposition fails to 
solve the separability problem is given by the state p= \ (P ei <8> P ei + P 'x ® P x ), 
where \x) = (|ei) + \c2))- The singular value decomposition of p is given by 

3 1 
p= -\u 1 ){u- l \ + -\u 2 ){u2\ (5.16) 

(f\ f K\ 

where = 2 _^_ and |u 2 ) = I f I . Clearly mat\u 2 ) ^ 0, thus we 

VSJ V-l) 

cannot conclude the separability of p. This illustrates the need to find the 
decomposition of p into rank-one partial isometries without the orthogonality 
relation. 



6 Conclusions 

In summary, we have reformulate the separability problem of quantum states as 
a descomposition into partial isometries of the density matrix transformed un- 
der the tilde transform. We have deduced several facts from this reformulation: 
(i) given that both the tilde transform and the singular value decomposition 

3 This follows from the preceding proposition and the fact that is is impossible to find four 
two-dimensional pairwise orthonormal positive matrices. 
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of an arbitrary matrix are computationally accesible tasks, we have provided 
a systematic sufficient test of separability both for bipartite and multipartite 
systems of any dimension, and (ii) we have also provided a complete solution 
independent of the Schmidt decomposition for the pure-state case in both bi- 
partite and multipartite systems. 

Nevertheless there is still a relevant open question: a criterion to discern 
whether a given matrix admits a decomposition into rank-one partial isometries 
fulfilling the conditions of theorem^or of theorem[21 This demands a systematic 
study of the cone of matrices whose extremals are of the form with 
\u) and \v) with structure |a;*} <g> \x). Although this remains to be done, the 
approach presented herein shows the advantage of being easily extended to the 
multipartite case. 
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